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Abstract 

We extend the quasilocal formalism of Brown and York to include electromagnetic and dilaton 
fields and also allow for spatial boundaries that are not orthogonal to the foliation of the spacetime. 
The extension allows us to study the quasilocal energy measured by observers who are moving 
around in a spacetime. We show that the quasilocal energy transforms with respect to boosts 
by Lorentz-type transformation laws. The resulting formalism can be used to study spacetimes 
containing electric or magnetic charge but not both, a restriction inherent in the formalism. The 
gauge dependence of the quasilocal energy is discussed. We use the thin shell formalism of Israel to 
reinterpret the quasilocal energy from an operational point of view and examine the implications 
for the recently proposed AdS/CFT inspired intrinsic reference terms. The distribution of energy 
around Reissner-Nordstrbm and naked black holes is investigated as measured by both static and 
infalling observers. We see that this proposed distribution matches a Newtonian intuition in the 
appropriate limit. Finally the study of naked black holes reveals an alternate characterization of 
this class of spacetimes in terms of the quasilocal energies. 



1 Introduction 



Gravitational thermodynamics continues to be one of the most active areas of research in gravi- 
tational physics. The interest began in the early 1970's with Bekenstein's recognition that if the 
temperature of a black hole is proportional to its surface gravity and its entropy is proportional to 
the surface area of its horizon, then the laws of black hole mechanics are laws of thermodynamics 
|TJ . These speculations were confirmed by Hawking's discovery that a black hole emits radiation as 
a perfect black body with temperature proportional to its surface gravity 0] and by calculations 
based on Gibbons and Hawking's proposal of the Euclidean path integral formulation of gravity^] 
which predicted that a black hole has an entropy equal to one quarter of its surface area. 

In the quest for a theory of quantum gravity, these laws of black hole thermodynamics are 
amongst the few solid results available. As such the investigation of their features and the search 
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for a statistical mechanics underlying them have been areas of tremendous activity and interest 
ever since. Indeed, the entropy/area relationship is viewed as an acid test which all aspiring 
quantum gravity theories must pass. In particular, in recent years there have been successful 
attempts to show that both string theory Q and canonical quantum gravity H correctly predict 
the entropy of black holes. 

On a more prosaic level however, one of the most popular current approaches to semi-classical 
gravitational thermodynamics is based on the Euclidean path integral formalism and a Hamiltonian 
analysis of the action functional that was proposed by Brown and York || They considered a 
spatially bounded region of spacetime and defined a stress energy tensor over the history of that 
system's boundary as the functional derivative of the action with respect to the three metric on that 
boundary. Foliating the spacetime they then defined quasilocal energy (QLE), angular momentum, 
and spatial stress densities as projections of the stress energy into the corresponding foliation of 
the boundary. Reference terms subtracted from these densities are defined by embedding that 
same boundary into a reference spacetime and recalculating the quasilocal quantities for that 
new embedding. Integrating the densities over one leaf of the foliation of the boundary, the full 
quasilocal quantities (rather than just tensor densities) are defined for that "instant" of time. 

With the energy, momentum, and pressure defined on the boundary one can do gravitational 
thermodynamics (see as examples @, || [| ) , study quantum mechanical black hole pair creation [ 10 
11, 12 1, and investigate the distribution of gravitational energy in a variety of spacetimes [13, Il4| - 



matter of interest in its own right. It has been shown that in the appropriate limits, the QLE agrees 
with many of the usual definitions of the total energy of the system. For example in asymptotically 
flat spacetime it has been shown to be equivalent to the Arnowitt-Deser-Misner (ADM) energy 
|| (in §) and the Trautman-Bondi-Sachs (TBS) p] energy (in |l^l). In asymptotically anti-de 



Sitter spacetimes it has been shown [g, [18| to be equivalent to the Abbot-Deser energy JT9[ . Very 



recently interest in this subject has been renewed with the AdS/CFT inspired redefinition of the 
reference terms with respect to intrinsic rather than extrinsic curvatures. The formalism can then 
be used to investigate the thermodynamics of a new group of (mainly AdS) spacetimes [20] that 
were inaccessible to the embedding reference term approach. 

Within the original quasilocal formalism however, there was an acknowledged incompleteness. 
Namely, it was assumed that the spatial boundary was always orthogonal to the spacetime fo- 
liation. This simplified some calculations and was true for all of the standard examples where 
one considers static spherically symmetric boundaries in static spherically symmetric spacetimes 
foliated according to the usual time coordinate. From a theoretical standpoint this assumption 
restricted the variations used to define the stress tensor (and calculate the equations of motion) 
to those that preserve the orthogonality of the foliation and the boundary. Practically, it meant 
that it was extremely difficult to calculate quasilocal quantities measured by non-static observers. 
As an example, the history of a spherical set of observers falling into a Schwarzschild black hole 
would not be orthogonal to the foliation defined by the usual time coordinate t. 

Only a few papers have considered the consequences of dropping the orthogonality restriction. 
Hayward [21| modified the gravitational action so that its variation is well defined for arbitrary 
variations of the metric. Lau partially dealt with the issue in [22] though his main concern was 
the definition of quasilocal quantities with respect to Ashtekar variables. Hawking and Hunter 
developed a non-orthogonal formulation which explicitly depended on the angle of intersection 
between the foliation and the boundary [23|. This dependence could only be removed by carefully 
foliating the reference spacetime to have the same intersection angles at the boundary, making the 
calculations quite cumbersome. 

Our own considerations on this issue were the subject of a recent paper [^4| . In contrast with the 
approach of Hawking and Hunter we focused on the foliation of the spatial boundary rather than 
that of the spacetime as a whole. With this approach the formalism was independent of the foliation 
of the rest of the space-time even without the inclusion of reference terms. Computationally it 
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was easier to apply and we also argued that our approach was the more natural one to take (see 
1 24] or section ||| of this paper for more details). In our approach it is relatively easy to calculate 
such quantities as the quasilocal energy measured by observers falling into a black hole. 

An interesting class of black holes, almost tailor-made to be investigated by our methods, has 
recently been discussed by Horowitz and Ross [25]. These so-called naked black holes are massive 
near-extreme static black hole solutions to the equations of string theory in the low energy limit. 
Observers who aren't moving with respect to such holes and who are just outside the event horizon 
measure only small curvature invariants. However, observers falling into the holes along geodesies 
see Planck scale curvatures and experience correspondingly massive crushing tidal forces as they 
approach the horizon. The holes were dubbed "naked" because Planck scale curvatures can be 
observed outside of their event horizons. 

In this paper we extend our previous work on non-orthogonal boundaries (reviewed in section 



3.l|) to include appropriate matter fields (section p.2|) so that we can calculate the QLE measured 



by these two sets of observers and compare it to the curvatures that they observe. The appropriate 
matter fields in this case are an electromagnetic field coupled to a dilaton field. As such along 
the way we will examine the gauge dependence of the quasilocal energy, see why this gauge 
dependence arises, and investigate exactly how it manifests itself. We shall see that the quasilocal 
energy naturally breaks up into a gauge independent geometric part and a gauge dependent part. 

The naked black holes that we wish to study are magnetically charged so we necessarily also 
confront issues related to electromagnetic duality (section |3.3j ) . We will show that the formalism 
that we have set up does exhibit this duality, but at the same time is not suited to dealing with 
dyonic black holes. In its current formulation, it can handle situations with electric or magnetic 
charge but not both simultaneously. 

Having dealt with these gauge-theoretic issues which are essentially unrelated to our non- 
orthogonal extension we then examine in some detail how quasilocal quantities transform with 
the motion of sets of observers who are measuring them (section [O]). We significantly improve 
our earlier treatment [24] and see that the quasilocal energy in particular transforms according 
to simple Lor entzian- type laws (though in this case there are two different velocities involved in 
those laws). 

Then, before turning to the examples, we examine the very close correspondence between the 
quasilocal formalism and the work by Israel explaining the behaviour of thin shells of matter in 
general relativity (section |3.5| ). This correspondence provides support for the quasilocal notion of 
energy and enables us to reinterpret it in an operational way. In the light of this reinterpretation 
we then examine recent AdS / CFT inspired work [^] on reference terms for the quasilocal energy. 

The first spacetimes that we investigate in the light of the preceding observations are Reissner- 
Nordstrom spacetimes (section [4.1[ ). We calculate the static and infalling quasilocal energies and 
the Hamiltonian for such spacetimes. We show how these quantities correspond to the energies 
that we would expect using a quasi-Newtonian intuition. We also see that infalling observers will 
measure arbitrarily large total energies as they fall into black holes that are arbitrarily close to 
being extreme. This contrasts with observers measuring the geometric part of the energy, who 
will find results only slightly larger in magnitude than those seen by their static counterparts. 

Having used the Reissner- Nordstrom spacetimes to orient ourselves and gain some intuition, 
in section 4.2 we finally apply the non-orthogonal formalism to our intended target - the naked 
black holes. There we find that the behaviour of quasilocal quantities calculated for these black 
holes is qualitatively the same as for near extreme Reissner-Nordstrom holes with one exception. 
For reasonable choices of the coupling constant between the electromagnetic and dilaton fields, 
static observers measuring the geometric energy will record large values (proportional to the mass 
of the hole) while those who are infalling will actually measure arbitrarily small energies as they 
cross the horizon. We discuss this rather surprising result in the final section of our paper. 

First though, as a preliminary to all of the above work, we establish a set of definitions and 
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examine the regions of spacetime that we shall be working with (section |2.1[ ) and then review the 
field equations of electromagnetism coupled to a dilaton field on a general relativistic background 
(section |2.2| ) . 



2 Definitions and Set-up 

Apart from some minor changes, most of the notation for this paper will be familiar to those who 
have read || and [ 24 1 . In this section we set up the notation and review some facts on matter fields 
that will be pertinent to the rest of the work. We begin by setting up the geometrical background 
in which we will be working. 



2.1 The geometry 

Let Ai be a four dimensional spacetime with metric tensor field g a p and in that spacetime define 
a smooth timelike vector field T a and a spacelike three dimensional hypersurface T,q. This field 
and surface are sufficient to (at least locally) define a notion of time over Ai; if Eo is taken to be 
an "instant" in time then past and future "instants" may be constructed by evolving So with the 
flow defined by the vector field T a . The resultant time foliation of Ai may be parameterized by a 
coordinate t such that T a d a t = 1. Then in the usual way we say that an "instant" £fi "happens" 
before £^ if ^i < t%- If we consider T a to guide the evolution of a set of observers in Ai then this 
is the observer defined notion of time. 

We may break up T a into its components perpendicular and parallel to the £t by defining a 
lapse function N and a shift vector field V a so that 

T a = Nu a + V a , (1) 

where u a is defined so that at each point in Ai it is the future pointing unit normal vector to the 
appropriate hypersurface £f The shift vector is constructed so that V a u a = 0. The lapse and 
shift then tell us how observers being swept along with the time flow T a move through space and 
time relative to the foliation. 

If we consider a set of observers forming a closed two surface f2o m then we may naturally 
define a four dimensional region M £ Ai. Let B be the timelike surface defined by evolving 0,q 
using the vector field T a . We use the foliation of M to induce a corresponding foliation of B; that 
is we define Slf = S^D-B. Assigning one side of to be an "inside" and the other an "outside", we 
may extend that notion over all of B. Finally, choosing £ti and £t2 as initial and final surfaces we 
may naturally define the region M £ Ai to be the region of spacetime inside B and that "happens" 
between times t± and £2- Figure [l| illustrates these concepts for a three dimensional Ai. 

Note that while a £t uniquely specifies an f^, the converse isn't true. Any number of £t 
foliations can be defined that are compatible with a given fit foliation. In fact despite the way 
that we have set up the foliations in this section, for most of this paper we will only be concerned 
with the foliation The foliation of the rest of the spacetime is irrelevant, basically because 
there are no observers in the interior of B to define it. The only observers reside on the boundary 
B. 

We define unit normal vector fields for the various hypersurfaces. Already we have defined 
u a as the future-pointing timelike unit normal vector field to the £t surfaces. Similarly, we may 
define u a as the future-pointing timelike unit normal vector field to the surfaces fit tangent to the 
hypersurface B. The spacelike outward-pointing unit normal vector field to B is defined as n a . 
Then, by construction u a n a = and T a n a = 0. We further define h a as the vector field defined 
on B such that h a is the unit normal vector to tangent to (Qj being viewed as a surface in 
St). By construction u a h a = 0. 
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We define the scalar field r/ = u a n a over B. If rj = everywhere, then the foliation surfaces 
are orthogonal to the boundary B (the case dealt with in refs. [l^] Q), and the vector fields with 
the tildes are equal to their counterparts without tildes. If 77 7^ 0, we express u a and n a in terms 
of u a and n a (or vice versa) as, 

n a = —n a — r\u a and u a = —u a — r]n a , (2) 
A A 
or, 

h a = \n a + V u a and u a = \u a + r]n a , (3) 
A A 

where A 2 = t-ts- 

Note too that on the surface B we may write, 

T a = Nu a + V a , (4) 

where we call N = XN the boundary lapse and V a = a^V^ the boundary shift. This is possible 
because we have assumed that T a n a = on B. 

Next consider the metrics induced on the hypersurfaces by the spacetime metric g a p. These may 
be written in terms of g a p and the normal vector fields. h a p = g a /3 + u a Uf3 is the metric induced on 
the E t surfaces, 7 a/ g = g a p — n anp is the metric induced on B, and o a p = h a p — n a np = j a f3 + u a U/3 
is the metric induced on Qf. By raising one index of these metrics we obtain projection operators 
into the corresponding surfaces. These have the expected properties: h a piiP = ■y a /3 n /3 = a a ^nP = 

a a p u^ = 0, and h a ph% = h a 7 , 7^7^ = 7 Q 7 , and a^a^ = a\ 

On choosing a coordinate system {x 1 ,^ 2 ,^ 3 } on the surface £0 we define h = det(h a p) (where 
in this case we take h a g as the coordinate representation of that metric tensor). We then map 
this coordinate system to each of the other Et surfaces by Lie-dragging with the vector field u a ; 

1 The definitions of u a and n a are consistent with ref. ^j, but interchanged with respect to those in ref. |]l~8| . 
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combining this set of coordinates on each surface with the time coordinate x° = t we have a 
coordinate system over all of M. We define g = det(g a p). Similarly, choosing a coordinate system 
on Qt we define a = det(a a g). Again, using the time flow to extend the coordinate system over 
all of B, we define 7 = det{^ a p). It is then not hard to show p3[ that 
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NVh and ^=7" = N yfa. (5) 



We also define the following extrinsic curvatures. Taking V Q as the covariant derivative on 
M. compatible with g a p, the extrinsic curvature of T* t in Ai is K a p = —fi' a h & pV^us = —\£ u haf3, 
where £ u is the Lie derivative in the direction u a . The extrinsic curvature of B in A4 is Q a p = 
— 7'a7 <5 flV 7 n5 while the extrinsic curvature of Of in T, t is k a p = —a^a^gVyhg. Contracting each 
of these with the appropriate metric we define K = h al3 K a p, = 7 a/3 Q/ 3, and k = & al3 k a p. 

Finally, we define the following intrinsic quantities over A4 and S^. On M, the Ricci tensor, 
Ricci scalar, and Einstein tensor are TZ a p, 1Z, and G a p respectively. e a g^s is the completely skew 
symmetric Levi-Cevita tensor. On St, D a is the covariant derivative compatible with h a /3, while 
R a p and R are respectively the intrinsic Ricci tensor and scalar. 

2.2 Fields on the spacetime 

We consider spacetimes containing a cosmological constant A, a massless scalar field <f> (the dila- 
ton), and a Maxwell field F a g. The field equations are: 



l -e a ^V p F l5 = 0, (6) 

V^e-^F^) = 0, (7) 

V a V a </> + ^ae' 2 ^ 'F a(3 F a P = 0, and (8) 

Gap + A gaf) - torTaf, = 0, (9) 

where a is the coupling constant between the scalar and Maxwell fields, and 

T at3 = i- ([V Q #V^] - l[V^)[V^]g a p + e- 2a *[F aj F^ - \g a pF l5 F~> 5 ]\ (10) 

is the stress-energy tensor associated with the matter. The first equation holds because the 
electromagnetism is described by a gauge field, while the last three may be derived from the 
action principle as we will do in section ||. However we first review some useful facts regarding 
these fields and their projection into three dimensional hypersurfaces. 

First define the dual *F a3 = \er 2a<i> e a ^ 5 F 7<5 of F aj3 . Then, we may respectively rewrite the 
above equations as 

V 3 {e 2a<t> *F al3 ) = 0, (11) 
-\e a ^ S Vp*F 1& = 0, (12) 

V a V Q </> - ^ae 2a * *F a3 *F af3 = 0, and (13) 
G a p + Kg ap - 8nT a3 = 0, (14) 

where this time we express the stress energy as 

T a3 = ^ (jy a <f>Wp<f>] - l^Wy^gap + e 2a ^F ai - ±g a/3 *F lS *F7*A . (15) 
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Thus, the equations of motion as a set are invariant under the duality transformation ((/> — > 
-<f>,F a p -> *F Q , J g). 

Second we define projections of the Maxwell Z7"(l) field into the spacelike hypersurfaces St 
considered above. If u a is the normal to that surface, then the (dilaton modified) electric and 
magnetic fields in that surface are E a = e~ 2a ^F a pu^ and B a = — '^e OL ^ <S vPF^s- Then, a simple 

calculation shows that E a = —\^ a p lS u^-kF^s and B a = —e 2a ^-kF a pu^ . Conversely we may write 
F a p and *F a p in terms of the electric and magnetic field three-vectors as 

F af s = e 2a *(u a E -upE a ) -e^B^ug and (16) 
*F a0 = - e - 2 ^{u a Bp-upB a )-t a /E lUs (17) 

Combining these two versions of F a p and *F a p with equations (Q) and ([ll]) and decomposing into 
the components parallel and perpendicular to the £j surfaces one recovers the (dilaton modified) 
three dimensional source-free Maxwell equations. In particular the components parallel to u a give 
us the divergence relations 

D a B a = and D a E a = 0. (18) 

Note that in terms of the electric and magnetic field, the duality transform F a p — ► *F a p becomes 
E a — ► —B a and B a — ► E a . 

Finally F a p is a 2-form field and by equation @ its exterior derivative is zero. Therefore (at 
least locally) there exists a one form field A a such that F a p = d a Ap — dpA a - This field is the 
electromagnetic potential. We decompose it into its components parallel and perpendicular to Y* t 
as a scalar potential $ = —A a u a , and three- vector potential A a = haAg, where $ is the Coulomb 
potential and A a is the vector potential. Then the electric and magnetic fields may be written as 

E a = -e" 2 ^ (JjD a (N$) + £ u A a ^j and (19) 
B a = -e a /u^D^As. (20) 

Similarly by equation (a) there locally exists a 1-form *A a such that -kF a p = d a -kAp — dp-kA a . We 
can decompose it into components = — * A a u a and -kA a = ha *Ab 0- Then in terms of these 
potentials the electric and magnetic fields may be written as 

E a = -e ocp lS u l3 D 1 *A 5 and (21) 

B a = e 2a * (^D a {N*$) + £ u *A a y (22) 

3 Quasilocal quantities from the action principle 

In this section we start from an action principle and define quasilocal energy, momentum, and 



stress tensor densities in the spirit of ||, 24 ] . We begin by reviewing the matter free case discussed 
in [24] and then include coupled dilaton and electromagnetic fields. We then see that the formalism 
that we have set up does not allow for magnetic charges and so use electromagnetic duality to 
modify it so that we can study magnetically charged spacetimes. Finally in the last part we see 
how the quasilocal quantities transform with the motion of the observers. 



2 There is an abuse of notation here. *A a , *A a , and *$ are in no sense duals of their unstarred counterparts. The * 
simply indicates their relationship to *F a p. 
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3.1 The matter free case 



Given M C M with the non-orthogonal boundaries described above an appropriate action for 
pure gravity (with a cosmological constant) is pi]]: 



I = — f d A Xyf^(K - 2A) + - / d 3 xVhK - - f dPxyf^Q 
2k Jm k is K Jb 

+ — / d 2 XA/asinh _1 (77) — /, 
K Jn 



(23) 



where J" s = and j n = — , and if we choose a system of units where c = G = 1, 

k = 8tt. In general / is a functional of the boundary metrics on dM. It will be discussed in more 



detail in section 3.4; for now we take it to be zero. 



As was shown in p4|] the variation of I with respect to the metric is, 

81 = f d 4 x^(G af3 + Ag a p)5g a(3 
Jm 



(24) 



f d 3 X (pfSKp) + ! d 2 X (p^ 



dt I d 2 x\fa 



N 



e5N -j a 5V^ - — 8^ ou aP 



In the above, P, 



aj3 



h 



s/h 

2, 



}(Kh<*P-K a P), P^ = isinh-S, e = \% j a = ±a£u 5 Vpn s , and 



s a/3 = \{k a P - [k 



2^V /3U a is the acceleration vector associated with u a . 



In the usual way if we solve 51 = while holding the boundary metrics constant, we recover 
the Einstein equations. Examining the initial and final hypersurfaces Si and £2 and their bound- 
aries fii and Pyf is the T, t hypersurface momentum conjugate to h a p while P is the Qf 
hypersurface momentum conjugate to \fo. Further, —\foi. is conjugate to the boundary lapse N, 
x/aja is conjugate to the boundary shift V a , and ^\fas a P is conjugate to the boundary metric 
(J a p. Following the Hamilton- Jacobi analysis of || we identify these three quantities as the sur- 
face energy, momentum, and stress densities. Note that they depend only on the foliation of the 
boundary and are indifferent to the foliation of the spacetime as a whole. 

I may also be decomposed with respect to the foliation to obtain pj] 



/ = f d A x (p aP £ T h a( 3 -NH- V a H a ) + [ dt [ d 2 xP^(£ T ^) 
Jm ^ 'J Jn t 



(25) 



dt / d 2 xyfa [Ne -V a j c 



Ti and TC a are the Hamiltonian and momentum constraints for the gravitational field on the £4 
surfaces and are defined as usual by 



Vh, 



Vh 



H = - — (G a p + Ag af3 )u a u l3 = -?-(R-2A + K 2 -K a pK a n=0, (26) 
and 

U a = —hfiGpy + Ag Pl W = — (OpRf - D a k) = 0. (27) 



They are zero for solutions to the Einstein equations. 

The quasilocal quantities are exactly those that would be obtained by the analysis of § if 
the foliation of B had been induced by a foliation £t of M that was perpendicular to B. 
Consequently we may think of them as being defined for a (local) orthogonal extension of the 



foliation of B into M. These are the natural quantities that a set of observers restricted to the 
surface B would measure. As was mentioned earlier, such observers know about the foliation of 
the boundary (it is defined by their notion of simultaneity) but being restricted to B they have 
no unique way of associating that foliation with a foliation of M as a whole. The natural foliation 
for them to work with is therefore the (local) orthogonal extension of B. In turn, this means that 
they will measure the quantities that we have seen arise from the action in a natural way. 

An observer-dependent Hamiltonian may also be defined. Recall that in elementary classical 
mechanics with one degree of freedom, the action I and Hamiltonian H are related by the equation 
I = J dt(pq — H), where q is the configuration variable for the system and p is the momentum. 
Extending this definition of the Hamiltonian to the system under consideration we obtain the 
Hamiltonian 

H= I d 3 x[NH + V a H a ]+ [ d 2 x^/a(Ne-V a j a ), (28) 
Js t Jn t 

for a spatial three surface bounded by Qf If T a is a Killing vector of the boundary metric "fa/3, 
then this Hamiltonian is a fixed charge. It is usually taken to be the mass contained within 
H, |J. Note that for solutions to the equations of motion the bulk constraint terms are zero and 
so the Hamiltonian depends only on the foliation of the boundary. 

Finally, before we move on to include matter fields, note that the above derivation assumes 
that within M the metric tensor is non-singular. If this is not the case, then the frequent uses of 
Stokes theorem in the derivation don't apply. Of course, black holes do have singularities in their 
centers and so problems arise if we wish to study them. For black holes, the full analysis only 
applies if B includes an inner boundary B' cutting off the singularity - in which case we must 
also consider the boundary terms on B' and are in fact considering the energy contained between 
the two bounding surfaces rather than in the hole itself. On the other hand if we consider a star 
there are no such singularities and it is reasonable to consider just the outer boundary B. Thus, 
in order to make a comparison between stars and black holes (which after all are described by 
the same solutions to the Einstein equations once you are beyond their defining surfaces) it is 
conventional to go beyond the preceding derivation and define quasilocal energies (and momenta) 
for black holes using just an outer boundary B. This is equivalent to the way one defines the 
electric charge of a point charge (as opposed to a charge distribution) in electromagnetism. 

There is an alternative but entirely equivalent way to look at this for spacetimes with an 
asymptotic region. Let us assume that the quasilocal energy E^ in the asymptotic limit is the 
total energy in the spacetime. As pointed out in the introduction, in this limit it does agree with 
most of the popular ways of defining total energy and so this is a reasonable assumption. Then for 
a given surface f?t in a leaf £( we may calculate the QLE contained between and the asymptotic 
surface at infinity. It is E^ — Eq. This region has no singularities to cause trouble and so the 
analysis may be executed rigorously. Next, the quasilocal quantities are additive so the energy 
contained within is — (E^ — Eq) = Eq which is the regular QLE. As noted the assumptions 
behind this approach are completely equivalent to those behind the one we discussed above, but it 
does give us a slightly different way of looking at things and highlighting what those assumptions 
really are. 

3.2 Including coupled electromagnetic and dilaton fields 

The above analysis may be extended in a straightforward matter to include gauge and dilaton 
fields. This extension (for the orthogonal case) was made quite generally in ref. || but for our 
purposes we just need to consider the coupled Maxwell and dilaton fields governed by the following 
action: 

lEMdii = 7^/ d A x^(K-2A-2(V a <P)(y a <P)-e- 2a <f > F Q(3 F a ( 3 ) (29) 
^ Jm 
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+- / cPxVhK-- [ d 3 xy/^/Q + - [ d 2 xyfc smb.' 1 (ri)- I. 
«is k Jb k J n 

For the rest of this section we shall assume that there exists a single (but not unique) vector 
potential A a properly defined over all of M such that F a g = d a Ap — dpA a (the converse of this 
assumption is that multiple vector potentials are required to generate F, each of which is defined 
over only part of M) . As we shall see in section |3.3| this is equivalent to assuming that there is no 
magnetic charge contained by any closed surface in M. 
Then, it is a simple calculation to show that 



6 (- 2v ^(V Q <A)(V» - y^e-^FapF * 1 ) (30) 
- V=sV a ([V a m + e-^F^SAp 



Further, 



Foil = V a V Q 0+iae- 2 ^F a/3 F^ (31) 

Fem = V aleT** F°0], and (32) 
2 f r _ „ r _ 1, 



la/3 



[V«0] [V^] - \ [V^] [V 7 0] g a ^j + e~ 2a *T™, (33) 



where T®Jf- = ■|(F Q , 7 F g 7 — \g a pF 1 sF~ i& ) is the standard electromagnetic stress energy tensor. 

The final term of equation ( |30| ) is a total divergence, and as such when we substitute it into the 
expression for SIemmi mav De moved out to the boundary using Stokes theorem (here we make 
use of the assumption that there is a single A a ). Then, using the vacuum result (|24|) we obtain 



SlEMdil = ^ 1 AV=0 {(G af 3 + Ag a p - 8irT aP )5g a ^ + VF Da 5<t> + ^ EM 5A^ (34) 
+ £ d 3 x [Pfdh a p + P^4> + P%8A a } + J d 2 x {P^(v^)} 

-JtoJ Q d 2 x^^(e + e^)5N-U a +%)5V a -^s a P5a a ^ 

+^JdtJ d 2 xy/^N {-£ n <t>6<t) + (npE^S® - e- 2a<t> u a n p e a/3 rS B^Ag} . 

In the above, Pf = ^{Kh ap - K aP ) and P^ = sinh ~ 1(??) as before while P^ = ^£ u 4> and 

P| = —^I±E a (recall that A a = h^Ap). They are the boundary momentum tensor densities 

conjugate to h a p, (j), A a , and yfa respectively, y/ae, y/crj a , and y/a^s a ^ are again the surface 
energy, angular momentum, and stress densities associated with pure Einstein gravity as defined 
in the previous section, while yfoe m and yfoj^ are the extra bits of surface energy and angular 
momentum density that must be added on to allow for the presence of the matter. They are 

£™ = ~np&k> and J^ 1 = -npE^A a . (35) 

<i> = —A a u a is the Coulomb potential and E a = e~ 2a ^F a pu^ is the electric field associated with 
u a , the timelike normal to in B. A a = a^A a is the projection of the vector potential into the 
Q t two boundaries and B a = — ^^^u^F^g. 
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Consider variations that fix h a p, <f>, and A a on the spacelike boundaries, and fix 7 a/ 3 (equiv- 
alently N, V a , cr a p), <f>, and jaAp (equivalently $ and aaAp) on the timelike boundaries. Then 
solving 5lEMd.il = we obtain the equations of motion (f?],||j9|) as promised in section 2J2 . Equiva- 
lently this particular action is only fully differentiate if we a priori fix all of the boundary metrics 
plus <ft and "faAp on the timelike boundary. 

It is not hard to see that fixing $ and Ap on the timelike boundaries is equivalent to fixing 
the component of B a perpendicular to B (that is B a n a ) and the components of E a parallel to 
Qt (that is apE 13 ). By Gauss's law the integral of B a n a over Qt is the magnetic charge contained 
by that surface. Thus the action is fully differentiable only if we restrict the parameter space of 
possible solutions to those with a specified magnetic charge which fits in nicely with our previous 
comment that the magnetic charge must be zero by our assumption that there exists a single vector 
potential generating the EM fields. In contrast, there is no restriction on the electric charge. We 
will further investigate this inequivalent treatment of the electric and magnetic charges in section 



3J. 

To decompose the action itself, we recalculate the gravitational constraints (this time including 
the matter fields) and obtain 

H m = - — (G a p + kg a p-MT aP )u a uP (36) 



H + — ([£P>][£> 7 <£] + {£ u ^f + e 2a ^E 2 + e-^B 2 



and 

n™ = —tfjGf h + h.g fh -%nTf h )u> (37) 
2\[h ( 



At the same time, starting from the definition of E a in terms of the potentials (equation 19) it is 
not hard to rewrite E a as 

E a = [D a [-N^> + V p A p ] - £ T A a + e^l^sVJ . (38) 

We may then combine these three results to decompose the action IemcLU as follows. First, after 
decomposing the purely gravitational terms as before we are left with the following as the bulk 
term integrand: 

P a/3 £ T h a0 - NH — V a H a - ^Z#V Q 0V> - ^-e- 2a(f, F a0 F al3 . (39) 

k 2k 

Bringing in (36) and (37) we may then rewrite this as 

P af3 £ T Kp - NH m - V a H^ + -^-(jV[i^] 2 + V a £ u 4>D a $) (40) 

K 

Next with ( |38| ) and the trivial £t4> = N £ u $ + V a D a (j) we may rewrite this entirely in terms of 
time derivatives of fields, constraint equations of those fields, and total divergences that may be 
removed to the boundaries. Then (59) becomes 



P aP £ T Kp + P^£ T ^ + P%£tA - NH m - V a H™ - T a A a Q (41) 
+ 2 -^D p (E^A a ). 
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Q = 2 -^D a E a is the constraint equation for the electric field with no sources (equation 18). Thus 
the final result is 

lEMdii = J dt J x d * x {pf£TKp + P^£ r + P^tAo} (42) 
+ j&j Q d 2 x [P^(£ T V^)} 
- J dt J d 3 x {NH m + V a H™ + T a A a Q} 

~ J dtj^ d 2 x^{N(e + e m )-V a (j a +^)}. 

Note that once again we have used Stokes theorem and therefore the assumption that there exists 
a single A a defined over all of M. 

Just as in section (|3.1|) we may define a Hamiltonian. Including the matter fields it is 



l EMdil 



H m =f d 3 x[NH m + V a H™ + T a A a Q]+ I d 2 x^ [n{e + e"») - V a (j a + %)} , (43) 
J St J fit 

for a spatial three surface St bounded by Qf Again however, the Hamiltonian is actually inde- 
pendent of that foliation Sj for solutions to the equations of motion. 

Note that even though the action lEMdii is g au g e invariant, this Hamiltonian does not necessar- 
ily inherit that invariance. The paths by which this gauge dependence can creep in are quite easily 
found but the effect is quite important so we will pause here to point them out in some detail. 
First, it is a simple matter to rewrite the action in terms of kinetic and Hamiltonian (potential) 
terms. Specifically, 

J dt |jf ^x(Pf£ T h a p + P^£T<f> + PX£TAa) (44) 

+ jf d 2 x(P^£ T ^) - tf m j - I. 

Thus while the time integrated sum of the Hamiltonian and kinetic terms must be gauge invariant, 
that invariance isn't necessarily inherited by the time integral of Hamiltonian itself unless part of 
the gauge freedom is used to ensure that £rA a = 0. If this is the case and H m is independent of 
the leaf of the foliation then H m will be independent of the remaining gauge freedom. That gauge 
and foliation are the natural ones to choose in stationary spacetimes. Nevertheless we should keep 
in mind that a (partial) gauge choice has been made. 

In the conventional usage of this work, there is an alternate route by which gauge dependence 
can also work its way into the Hamiltonian. Namely, as was discussed at the end of section EO , 
if we consider spacetimes containing singularities then the above analysis only strictly applies if 
we include an inner boundary to cut out that singularity. For non-singular spacetimes, there is 
no need to include this second boundary however, and so in the interest of comparing singular to 
non-singular spacetimes, it is conventional to work with the outer boundary only. Without that 



inner boundary however, the gauge dependence returns. In section ^4.1.4! we see how this shows up 
in a Reissner-Nordstrom spacetime. There we will also see that the remaining gauge dependence 
resulting from neglecting the inner boundary amounts to little more than a choice of where to put 
the zero of the electromagnetic energy. 

Thus, we see that the gauge independence of the action doesn't necessarily ensure the gauge 
invariance of the Hamiltonian. Specializing to the case where the lapse N = 1 and shift V a = 
on the boundary, and leaving aside the issue of how such a choice affects the relative foliations of 
the inner and outer boundaries, we note that the QLE will not in general be gauge independent 
either. 
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3.3 Electromagnetic duality 



In the previous section we started out by making the assumption that there was a single vector 
potential A a defined over all of M. This assumption immediately means that there is no magnetic 
charge in M (or contained by any surface residing in M). To see this let tlx be any closed spatial 
two surface in M with normals u a and n a . Then, the magnetic charge contained within fix is 
J Qx d 2 x^n a B a . By equation (£(]), n a B a = -rPvPe^ 6 D^As = d* \-n a vPe a ^ %) where d* 
is the covariant derivative in the surface fix- But this is a total derivative and so integrated over 
a closed surface it is zero 3 . Thus there is no magnetic charge contained by any surface in M. 
Thinking back to the discussion following equation (^) this agrees with the observation that the 
magnetic charge must be fixed under variations. Indeed it must, and that fixed value is zero. 
There is no room for magnetic charge within the set of assumptions that we have made. 

Keep in mind that this is a stronger statement than just the local statement that F = dA =>■ 
dF = d(dA) = =>■ D a B a = 0. When working with a vector potential, the manifestation of 
magnetic charge in the potential is global and topological (resulting from a twist in the U(l) 
gauge bundle) rather than local as is the case for electric charge. If we assume that there is a 
single potential covering M then the U(l) gauge bundle is trivial by definition and so there is no 
magnetic/topological charge. Even more strongly, as noted earlier, no surface contained in M can 
contain magnetic charge. This means, for example, that if M is the region contained within two 
concentric spheres (multiplied by a time interval), then not only is there no charge in M but also 
there is no charge in the region inside the inner sphere. 

In fact, projecting into spatial slices E| of M, de Rhams theorem (see for example p(|) tells 
us that a single vector potential is defined over all of if and only if there is no magnetic 
charge contained within any two surface fix C Sj. Thus if we wish to allow magnetic charge in 
our spacetimes we must also break M into at least two regions each of which has its own vector 
potential. Then, the frequent uses of Stokes theorem in the derivation will remove total divergences 
to the boundaries of those regions rather than just the boundary of M itself. By definition some 
of those region boundaries will actually be interior to M and so observers inhabiting dM will not 
be in a position to measure all of the boundary terms and therefore will not be able to fully assess 
what is happening in the interior of M. 

We can however study magnetically charged spacetimes if we use a duality rotation to make 
the magnetic charge local and analytic (while the electric charge becomes global and topological), 
in which case the action IemcLU becomes 

*I E Mdil = ^ / d 4 x^(^-2A-2(V^)(V»-e 2a nF^*F Q/3 ) (45) 

+- / d 3 xVhK-- [ d 3 x^/^j& + - / d 2 x^asinh- 1 (ri)-L 
tJs k Jb K Jn 

Note that since F a pF a P = — * F a p * F a " this action is not equal to Iem<M- Nevertheless this is 
the correct action to use if we wish to obtain the usual equations of motion. To wit, assume there 
exists a single vector potential *A a that generates *-F a /? over all of M. Then, the variation of the 
action is 

Si+lEMdil) 

= J- / d 4 x^|(G Q/3 + A^-8vrr^)^ + 4*^ M ^ + 4*^ M 5(^)) (46) 

+ J d 3 X [pfShaf, + P^H + P^K)} + J d 2 X {P^5(y/j)} 

3 In the more efficient differential forms notation, in the spatial slice orthogonal to u a , A is a one form and B — dA is 
a two form. Consequently the magnetic charge contained within fix is J n B = J n dA — since fix is closed. 
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+ 



JdtJ^ d 2 x^ I (e + *i™)SN - (j a + *%)5V a - js a ^5a aP j 
^JdtJ d 2 x^N {-£ n <j>8<f> - (npB p )6(*ib) - e 2a< ^ a n /3 e Q/3 7<5 'E 7 6(*A S )} 



-I- / dt I d 2 x\P, 



where most of the quantities retain their meanings from (|34|) but the starred quantities have 
undergone the usual transformation. Ergo, *Fdu = —V a V a 4> + \ae 2a< ^ * * F a @ , *T^ M = 

V Q (e 2 ^ * F^), - ^A^, P£ = ^B a , *e m = -fn^ * $, and ★JJ = "K^) * A*- 
If the boundary terms are all zero we obtain the usual equations of motion. The quantities that 
must be held constant in order to get those boundary terms equal to zero are unchanged for the 
pure gravitational/geometric terms, but the electromagnetic terms have reversed so now E a n a 
(and consequently the electric charge) and (TpB^ should be fixed. Of course by our assumption 
about the vector potential there is no electric charge and so, as for the magnetic charge in the 
previous example, it is naturally fixed. 
The decomposition of this action is 

* lEMdii = J dt d 3 x {pf£ T Kp + P^£ T <\> + P? A £ T (*A a )} (47) 

- j dt j d 3 x {NH m + V a H™ + (T a * Aa) * Q} 
-JdtJ^ d 2 x^ [N{e + - V a (j a + *#)} . 

TL m and 7i™ are invariant under the duality transformation, but *Q = —^^ L B'a_B a . 

Thus we have well defined formalisms with which we can study spacetimes containing either 
only electric or only magnetic charges. We do not, however, have a formalism that easily handles 
dyonic spacetimes. Admittedly we could make a duality rotation of the action to study a spacetime 
with a particular dyonic charge but that would still not allow us to consider spacetimes contain- 
ing multiple dyons with varying ratios of electric and magnetic charges. Furthermore, there is 
something fundamentally unsatisfying about having the form of the action depend on the charges 
contained in the spacetime. As it stands we don't have a solution for this problem and so will not 
consider dyonic spacetimes in this work. 

Finally, we recall that the relationship between electric and magnetic charge and an F a pF a ^ 
type action has been considered before in studies of the production of charged black hole pairs 



1 27, ^g, 11, 12]. In these cases one computes the probability for an empty spacetime with a source 
of excess energy, such as a cosmological constant, to quantum tunnel into a spacetime containing a 
pair of black holes. This process is calculated within the path integral formalism where the actions 
of instanton solutions to the equations of motion are taken as the lowest order approximation of 
the full Euclidean path integral. Those instantons (and indeed the set of all paths integrated over 
in evaluating the path integral) are required to satisfy certain boundary conditions. First, they are 
required to interpolate between the initial and final spacetimes matching onto these along spatial 
hyper surf aces, and second, boundary conditions are enforced to keep appropriate thermodynamic 
(quasilocal) quantities constant. The action that one uses to calculate the path integral must be 
chosen so that its variation fixes those same quantities. For the creation of a pair of magnetically 
charged nonrotating holes the magnetic charge must be fixed while in the creation of a pair of 
electrically charged nonrotating charged holes the electric charge must be fixed. 

Thus, in those papers IemcHI is used to study the creation of magnetic black holes, while 
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IemcLU + AlEMdil 18 used to study electric black holes where 

AI EMda = -- [ d 3 xVhe- 2a *u a F a PAp + - [ d 3 x^e- 2a ^n a F^A p . (48) 

In line with this paper we have included a non-zero dilaton. Then the variation of Iem<m + AIemcLU 
fixes the electric rather than the magnetic charge - the variation of AIemcHI switches the boundary 
terms to those of I+EMdii- This isn't really all that surprising of course since for solutions to the 
(in this case dilaton- modified) Maxwell equations, 

AlEMdil = - [ d A x^e- 2 ^F a pF a P (49) 

K JM 

= J_ f d A x^-g(e^F a(3 F^ -e 2a ^F a ^F^), 

^ JM 

and so for solutions lEMdil + AI E Mdil = I+EMdil- 

In light of these observations we can clarify what is happening for various choices of the action 
in pair creation processes. Our approach is commensurate with that taken in the electric case: 
the correct action to start with is lEMdil- However, as noted above, variation of J-EMdil implies 
that the electric charge is not fixed, and so a boundary term AlEMdU must be added to correct 
this problem. This term is well defined since for a purely electric black hole there is a single A a 
covering all of M. In a similar way a boundary term would have to be added to fix the angular 
momentum of the hole if we were allowing for rotation. 

By the same reasoning, the correct action to start with for a magnetic black hole is I+EMdU- 
To fix the magnetic charge one would then have to add on a boundary term AI+EMdih defined 
analogously to AlEMdU- This is equivalent to using lEMdil m considering pair production because 
for solutions to the equations of motion I+EMdil + AI+EMdU = J-EMdil- However we note that for 
J-EMdil the derivation of the QLE from the action only properly goes through if there is a single 
A a defined over all of M, in which case there can be no magnetic charge. 



3.4 Transformation properties and reference terms 

In this section we consider how the quasilocal quantities defined above transform with respect 
to boosts of the observers. We will also see how the choice of a reference term / changes these 
transformation properties. We will only explicitly develop the transformation laws for the F a p 
formulation. The extension to the dual formulation is trivial. 

Consider a set of observers being evolved by the vector field T a and with an instantaneous 
configuration f^. To avoid unproductively untidy notation we will take B to be orthogonal to the 
foliation T, t . This does not reduce the generality of what follows since we have already seen that 
the tilded quantities defined with respect to the foliation of B are equivalent to untilded quantities 
defined with respect to a local foliation orthogonal to B. As before the foliation is evolved by a 
vector field T a = Nu a + V a . We next consider a second set of observers boosted with respect to 
the first set. They have the same instantaneous configuration Q t but instead are being evolved by 
a vector field T* a which may be written in terms of a lapse N*, shift V* a , and timelike unit normal 
u* a to Qf The spacelike normal for this set of observers is then n* a - defined by the requirement 
that n* a u* a = and crfn^ = 0. Recycling earlier notation, we define 77 = u* a n a and A = —j==. 

Then 

u *a = la + &nd n *a = 1 a + ( 5Q ) 

A A 

Next, consider a set of observers who are static with respect to the foliation and who measure 
time as their proper time (i.e. they are being evolved by the vector field T a = u a ). Define 
an orthonormal frame for each of these observers {u a , n a , ef, } where ef^ef €= T£lf Then 
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this set of observers sees the T* a observers as having velocity components v\- = — = 

vi = - r*^ e J° = j^V* a e la , and v 2 = = j^V* a e la in the directions n a , ef, and e£ 

respectively. 

It will also be convenient to define the extrinsic curvature of Qt in B as kj^n = —aZc^V^ug, 
which is contracted to k\ = cr^k^. The rate of change of n a in the direction it points is 
a} = n^V pn a . The projection of the electromagnetic vector potential is <J?J = —A a n a . The choice 
of the | superscript is meant suggest an interchange of u a and n a in these quantities (as compared 
to the same expression without the superscript). We then define versions of the quasilocal densities 
with u a and n a interchanged: 



-k\ (51) 

K 

(52) 

K 

A = -^ a n 5 V^u 5 , (53) 

K 

fA = l( n P E p)A a , and (54) 

K 

s l aP = l(kl p -[kl-u^a\]a a p). (55) 



The | notation has been slightly abused in the matter terms as the u a and n a were not interchanged 
in the npE@ terms. Note too that j\ = —j a and j™ = j™- 



Some of these quantities were first used in |22j in the context of defining quantities that are 
invariant with respect to boosts. The simplest example of such an invariant, modified for the 
electromagnetic field included here, is (e+e m ) 2 — (e^+e m ^) 2 . This is analogous to m 2 c 2 = E 2 —p 2 c 2 
which is an invariant for a particle with energy E and momentum p in special relativity. This 
suggests that we view el + e m I as a momentum flux through the surface Q,f Support for this 
interpretation comes when we note that 

^ E t = -^L^P £ u(T = _}i£ u ^. (56) 

That is, el is zero if and only if the observers don't see the area of the surface they inhabit 
to be changing. If the area does change then they see a momentum flux through the surface 
which is something we intuitively associate with motion. Note that this means that a sphere of 
observers moving at constant radial speed in flat space will measure a momentum flux so the 
notion isn't entirely in accord with intuition. Of course without reference terms such observers 
will also measure a non-zero quasilocal energy so this is not entirely unexpected. 

By a series of straightforward calculations we obtain expressions for the quasilocal quanti- 
ties seen by the T* a observers in terms of the quantities measured by the T a observers. These 
transformation laws are 

£ * +£ m* _ _I f a ^ Van * f3 + 2 A a u* a {n*PE*p) 

K \ 

= i(e + e m ) + r?(el+e m l), (57) 
A 

fc + jT = -Uu^V p n; + 2(n^E* p )A a 

K \ 

= tia+jZ)-^<£VpV, and (58) 

K 
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1 t A 

= T-Sa/3 + V 3 a/3 + -°"a/?u* 7 V 7 r/. (59) 

In certain cases, these laws simplify into very familiar forms. Assume that e* = e m ^ = and take 

T* a = u* a (that is the T* a observers move only in the direction n a perpendicular to Qt and use 
the proper time as their measure of time). Then, 

e* + e m * = j(e + E m ), (60) 

fa+Ja* = fa + Ja*, and (61) 

s *a0 = l s aP H CT a pU* J V 7 (^yv\- ) , (62) 



where 7 = , = = yl + rj 2 = j . So, in this case the energy density transforms as we would intu- 

itively expect from special relativity. The angular momentum density is invariant as we would also 
expect since the only motion is perpendicular to its direction. The stress tensor has a somewhat 
more complicated transformation law that is dependent on the perpendicular component of the 
acceleration of the second set of observers. Breaking it up into its shear rj a p = s\ a m and pressure 
p = cr al3 s a f3 components we obtain a little simplification. Namely, rf a ^ = jrj a p and no longer has 
an acceleration dependence. However, p* = \p + ^"U* 7 V 7 ?7 and the dependence remains. 

The transformation laws change if I / 0. Physically the reasons for this are as follows. Recall 



24] that the I 7^ may be viewed as a choice of the zeros of the quasilocal quantities. The usual 



way to set these zeros is to choose a reference spacetime (M, Qq) in which one demands that the 
quasilocal quantities all measure zero. Then we define 1_ as follows. First, one embeds ($lt 
into (M_,g a n) and defines a vector field T a such that: 1) T a T a = T a T a , 2) £Tcr a /3 = £T& a f3 (in 
the sense that their projections into Q t are equal), and 3) the projections of T a into Q t are also 
equal (that is the boundary lapses are equal). This will not necessarily be possible for all choices 
of {Q, t ,a a p) and T a , but will be possible in a wide range of cases. Note that satisfying these 
embedding conditions amounts to embedding (Q t , a a p) in (M_,g ) and requiring that T a and T a 

evolve (Qfi a ap) m the same way as in M. It is also equivalent to a local (in the coordinate time 
sense) embedding of B in M . 

Then, a typical definition for J is 

1= I dt I d 2 x^{Ni-V a j a ), (63) 
j j 

where i and j a are defined in the same way as usual except that this time they are evaluated for 
the surface Q,t embedded in the reference spacetime. Thus, as noted, the choice of I effectively 
sets the zeros for the quasilocal quantities. As required by the formalism, J is defined entirely 
with respect to quantities defined over B. 

Now, given T a and T a observers watching T* a and T_* a observers, in general r/ = u* Q n a will 
not be equal to r/ = u* a n a . Physically this means that in order for (f^, a a p) to evolve in the same 
way in the two spacetimes, that surface will have to "move" at different speeds in each of the two. 
Then the transformation law for the quasilocal energy density with reference terms becomes 

e* + e m * - e* = i(e + e m ) + r](e^ + e m I) - + vA . (64) 

Thus, even if the T a observers are static and the T* a observers move perpendicularly to as 
discussed in the earlier example, the simple Lorentz transformation law will not hold, as the 
basic and reference components of the QLE will transform according to different perpendicular 
velocities. This situation is reminiscent of the relationship between spacetimes which differ from 
one another by conformal transformations; in general the quasilocal mass for a system of finite 
size will not be conformally invariant due to the lack of conformal invariance of the background 
action functional which specifies a reference background spacetime p9[. 
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3.5 Thin shells - an operational definition of QLE 



In this section we examine in some detail a correspondence between the quasilocal formalism and 
the mathematics describing thin shells in general relativity which was developed by Israel in 1966 
1 30 ] - This was noted in passing in || but here we shall examine it in more detail and use it to 
reinterpret the quasilocal energy from an operational point of view. We begin by reviewing the 
purely gravitational thin shell work, then add in matter fields, and finally consider the implications 
of this correspondence for the recently proposed AdS/CFT inspired reference terms. 



3.5.1 Pure gravitational fields 

Israel considered the conditions that two spacetimes, each of which has a boundary, must satisfy 
so that they may be joined along those boundaries and yet still satisfy Einstein's equations. He 
showed that as an absolute minimum the spacetimes must induce the same metric on the common 
boundary hypersurface. Further the Einstein equations will only be satisfied at the boundary if its 
extrinsic curvature in each of the two spacetimes is the same. If those curvatures are not the same 
then a singularity exists in the (joined) spacetime at the hypersurface. However that singularity is 
sufficiently mild that it may be accounted for by a stress energy tensor defined on that boundary. 
The change in curvature may then be interpreted as a manifestation of a thin shell of matter. 

Modifying Israel's notation and sign conventions to be compatible with those used in this paper 
that stress tensor is defined as follows. Consider a spacetime Ai divided into two regions V + and 
V~ by a timelike hypersurface B. Let the metric on V + be g^n and the metric on V~ be g~g, and 
assume that they induce the same metric j a p on B. Further, let n +a and n~ a be the spacelike 
unit normals of B on each of its sides and define and 0"^ to be the extrinsic curvature of 
B in V + and V~ respectively. Then, Einstein's equation will only be satisfied if a thin shell of 
matter is present at B with stress-energy tensor S a p = ^ {(©^/j ~~ @ + 7a/?) — (®a/3 ~ © _ 7a/3)}- 
Note that this is the stress-energy tensor in the surface. To write it as a four dimensional stress 
energy tensor we must add in an appropriate Dirac delta function. 

Now let fit be a foliation of B generated by a timelike vector field T a = Nu a + V a (which 
as usual lies entirely in the tangent space to B). Then observers who are static with respect to 
the foliation will observe the thin shell to have the following the energy, momentum, and stress 
densities: 

(65) 

a^^V-fUj | , and (66) 
n +s as)a af3 ) - (k~ p - (k~~ - n' s d 5 )a a(3 )^ , (67) 

where k^n = —aZapVyrijf and = a*' 3 ^ are, as they were in previous sections, the extrinsic 
curvature of the surface fit in a (local) foliation of Ai perpendicular to B. a a retains its earlier 
meaning. 

The correspondence between the quasilocal and thin shell formalisms is now obvious. Consider 
the surface (B,^ a p) embedded in a spacetime (Ai,g a p) and a reference spacetime (A4, g „). Fur- 
ther let (Ai, g a p) be isomorphic to (V + , g^n) (or more properly the portion of (Ai,g a p) to one side 
of B is isomorphic to (V + , g^)), and in the same sense let (Ai, g„) be isomorphic to (V~,g^g). 
Then for observers living on B and defining their notion of simultaneity according to the foliation 
fit, 

£ = e-e, (68) 



£ = S^u 13 
Sap = S-iSvZvp 
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Jot = jot - j a , and (69) 
S a p = S aj 3 — s a/3 , (70) 

where s^p is denned in the obvious way. 

We may interpret this mathematical identity of the formalisms in a couple of ways. First, 
as in H we may note that the quasilocal work formalism provides an alternate derivation of the 
thin shell junction conditions and stress energy tensor. Namely we may consider two quasilocal 
surfaces on either side of the shell and consider the limit as the two go to the shell. In that case 
any reference terms will match and cancel and we will be left with the stress energy tensor defined 
above. This derivation is quite different from the one used by Israel. 

From a slightly different perspective we may view the thin shell work as providing an operational 
definition of the quasilocal energy. Given a reference spacetime M_ which is assumed to have energy 
zero, then the QLE contained within a two surface Of of a spacetime A4 can be defined as the 
energy of a shell of matter 4 in A4 that has the same intrinsic geometry as Of (including the rate 
of change of those properties - see the previous section) and a stress energy tensor defined so that 
the spacetime outside of O t is identical to that outside of Of in -M. while inside it remains Ai . 
In fact, the QLE with the embedding reference terms that we have considered is defined if and 
only if the fields outside of Of can be replicated by a shell of stress energy with the same intrinsic 
geometry embedded in Jvi . Provided that we can embed Of in the reference spacetime A4, by the 
construction that we have considered in this section we can define the relevant stress energy for a 
shell in M. 

This operational interpretation of the QLE serves to highlight certain aspects of its definition. 
First, in the spirit of the Gauss law from electromagnetism the quasilocal approach assumes that 
a bulk property of a finite region of spacetime is fully reflected in the values of fields measured 
on the surface of that region. Just as the Gauss law of electromagnetism yields a measure of the 
electric charge contained within a surface from the components of the electric field perpendicular 
to it, here we attempt to define the energy contained in a region from properties of its boundaries. 
Therefore two different configurations of fields and matter that produce a given set of fields at a 
given surface will - by the very nature of the approach - be defined to have the same energy. This 
correspondence between the quasilocal and thin shell approaches strongly support this assumption. 

Related to this is the fact that the quasilocal energy as defined by Brown and York (and 
extended here), is additive. That is, the total energy contained within regions O a and 0& is equal 
to the energy contained within O a U05. This is a very natural requirement for energy to meet and 
combined with a definition of the total amount of energy in a spacetime (such as the ADM energy) 
is basically equivalent to the "Gauss law" assumption that we noted in the previous paragraph. 
If we know the total energy in a spacetime and also can measure the amount of energy outside of 
a surface Of then the additivity tells us how much energy is inside Of and also requires that any 
configuration of matter and fields that produces the correct geometry at the surface Of must have 
the same energy. 

3.5.2 Including matter fields 

So far we have only considered the purely gravitational case. With the inclusion of non-zero 
electromagnetic fields in the outer spacetime we must also include charge and current densities in 
the thin shell if the inner (reference) spacetime doesn't have such fields. These charge and current 
densities are defined to account for discontinuities in the electromagnetic field just as the stress 
tensor is defined to account for discontinuities in the gravitational field/geometry of spacetime. 
Their calculation is an exercise from undergraduate electromagnetism Specifically if there is 
no EM field within the shell then for the foliation Sf the electric charge density on the shell Of is 
ln a E a . 
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Given an EM potential A a , observers on the surface of the shell who evolve via the vector 
field T a will define a Coulomb potential —T a A a = — V a A a . In the usual way we then 
define the energy of the charge density in the field as the charge times the potential. That is, 
^n a E a (—T a A a ) = Ne m — V a j^. As usual this component of the energy is gauge dependent. 

The surface charges and currents do not change the definition of the surface stress energy 
tensor which was defined entirely by the Einstein equations. As such they also don't change the 
definitions of £, J a and S a p. Therefore, if we include the stress energy with the energy of the 
shell in the gauge field, the total energy density in a thin shell evolving by the vector field T a is 
N(e + e m ) — V a (j a This of course is exactly the same as the QLE of the region of space on 

and inside of the shell as measured by a set of observers being evolved by the same vector field, 
and so we see that the correspondence between thin shell and quasilocal energies remains. 

The preceding reasoning may be repeated for the dilaton. The dilaton charge in a given volume 
is given by the integral of n Q V Q </> over the surface enclosing that volume. For black hole solutions, 
the value of the dilaton charge is constrained by demanding the spacetime has no singularities on 
or outside of the outermost horizon p^] . In the thin shell case, n a V a (j) yields the dilaton charge 
density on the shell f2 t if there is a constant dilaton field inside. 

3.5.3 The AdS/CFT inspired reference terms 

Recently there has been quite a bit of interest (references |2(|) in defining the reference terms 
with respect to intrinsic quantities of the boundaries rather than their extrinsic curvature when 
embedded in a reference spacetime. These new reference terms have been inspired by the AdS / CFT 
correspondence and there is much to say for them. In the first place apart from the asymptotically 
flat cases it isn't really obvious what spacetime should be used as a reference. In the second, 
even after that choice has been made it is in general either computationally difficult or actually 
impossible to embed a given quasilocal surface in a given reference spacetime (as an example see 
[i~3H where Martinez tries to embed the an r = constant t = constant surface from Kerr space into 
flat space). Intrinsic reference terms remove this problem and as such are of great interest. 

Unfortunately however, the correspondence between the thin shells and QLE does not (in 
general) hold for the intrinsic reference terms proposed so far. For example consider the flat space 
reference used in some of these papers. There e = ^v 7 ^!, where is the Ricci scalar for 
fit- Recalling some elementary differential geometry, the extrinsic curvature of a two surface in 
flat three space is k = + while the intrinsic curvature is = x lx 2 ' wnere -^1 anc ^ -^2 
are the principal radii of curvature for the surface. Then, k > 2V i?( 2 ) and the equality only 
holds if X\ = X2 ■ That is the two reference energies are only exactly equal (and so the thin shell 
correspondence only holds) when fit is a sphere. However, as Lau pointed out j2(| the equality will 
also hold in the limit of an arbitrarily large fit in which case the difference between the two terms 
goes to zero. Thus, in the cases where we can compare to thin shell calculations the quasilocal 
energy with this intrinsic reference term only matches for spherical shells or in the asymptotic 
limit. 

Perhaps even more seriously the above observation implies that the energy contained by a finite 
ellipsoidal surface is non-zero in flat space while the energy contained by a similar sphere is zero. 
Given these observations it is probably best to follow Lau and view the these new reference terms 
as a short-cut for calculating the embedding reference terms in the asymptotic limit rather than 
as a replacement for them. 

4 Examples 

We now apply the foregoing work to a couple of examples. In the first we investigate the quasilocal 
energy distribution in a Reissner-Nordstrom (RN) spacetime. In the course of this investigation 
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we use the thin shell analogy to show that the quasilocal energy reduces to the intuitive New- 
tonian limit. For the Hamiltonian we also explicitly illustrate our earlier comments on gauge 
(in) dependence. Finally, in the last example we will come to grips with naked black holes. 



4.1 Reissner-Nordstrom spacetimes 

In this section we study the distribution of energy in a RN spacetime. The dyonic solution has 
metric 

dr 2 

ds 2 = -F{r)dt 2 + ^-r + r 2 {d6 2 + sin 2 6dip 2 ), (71) 

where F(r) = 1 - ^ + =o£a. m is the mass, and Eq and Go are respectively the electric and 
magnetic charges of the hole. The accompanying electromagnetic field is described by 

F = - -^dt A dr + G sin OdO A dip. (72) 
A local vector potential generating this field is 

A = -— dt -G cos edp + dx, (73) 
r 

where x = x{t^ r ^,<p) is any function over M. Note that the potential given above with x = 



is not defined over all of M since dip is not defined at = 0, ir. In section 3.3 we saw that our 
Lagrangian formalism as constituted is not suitable for discussing dyonic spacetimes. As such we 
focus on electric black holes in the following subsections. The results for magnetic black holes are 
identical if we switch Eq and Go in the following and add in the appropriate *'s. 

Setting Go = we calculate the quasilocal quantities measured by static, spherically symmetric 
observers. Specifically, we consider a surface of observers f^o defined as the intersection of the 
t = to surfaces r = ro surfaces (where to and tq are constants). They are evolved by the vector 
field T a = N(r)u a . N(r) is the lapse function while the shift V a = 0. For any choice of 
N(r) the observers will be static in the sense that they don't observe any changes in the spatial 
metric; the lapse just determines how they choose to measure their time. In particular, choosing 
N(r) = y/F{r) the observers measure time according to the coordinate t, while choosing N(r) = 1 
the observers measure time in the "natural" way (that is T a T a = — 1). 

Then, u a = -^=5f , n a = y/F(r)d?, and a series of straightforward calculations yields 

e m = _L ^(^o-r^) and 
4vrr 2 



1 

Airr 



(76) 
(77) 



We have made the substitution n = 8tt and so henceforth work in geometric units where G = c = 1. 
The angular momentum terms are all zero as we would expect for these non-rotating spacetimes. 



4.1.1 Static geometric energy 

Here we calculate the quasilocal energy associated with the density e. We label it the geometric 
energy since it depends only on the extrinsic curvatures. Then, not including the reference term 
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In the large r limit this becomes Eceo = —r + m + ^{m 2 — Eq). The e reference term is 

E= ( d 2 x^e = -r, (79) 

so 

E G eo ~ E = r - ^ 'r 2 - 2mr + E 2 ^m+—(m 2 - E 2 ), (80) 

in the large r limit. 

This is what we would expect from an application of Newtonian intuition to the (equivalent) 
thin shell situation. For this viewpoint consider how much energy it would take to construct 
a shell of radius r with mass m and charge Eq from material residing "out at infinity". Using 
Newton's and Coulomb's laws, it is straightforward to show that —^{rn 2 — E 2 ) units of work are 
required to assemble the shell. It is then very natural to say that this energy is "stored" in the 
field, outside radius r. Now, an observer sitting far from the hole and measuring how a neutral 
particle is accelerated by the gravitational field would say that the total mass contained in the 
spacetime is m. Then assuming that the energy is additive, the energy contained on and/or inside 
the shell with radius r is 

(total energy) — (energy in fields outside of r) = m + — (m 2 — Eq), (81) 

as we calculated above. An alternative way of looking at this is to say that m units of energy 
were used to create the mass m at infinity and then |(m 2 — E 2 ) units of energy were stored in the 
fields outside of the shell. Assuming that conservation of energy holds once the matter has been 
created then we again obtain the above result once the shell has been constructed. In either case, 
our calculation reduces properly in the Newtonian limit. This limiting case was first considered 
in the original Brown and York paper ||. 

Returning to the full expression we note that Eceo — E monotonically decreases as r increases, 
reaching a minimum of m at infinity. Thus the energy contained in the fields is negative, which is 
what we would expect for a binding energy. 

Finally note that for an extreme black hole where \Eq\ = m, EQ eo — E_ = m is a constant. From 
the Newtonian shell point of view this makes sense. If we consider the construction of a shell with 
charge \Eq\ = m and mass m out of particles which also have equal mass and charge, then equal 
but opposite electric and gravitational forces would act on the particles during the construction. 
Thus, no work must be done to build the shell and so no energy is stored on the fields. Alternately 
equal amounts of positive and negative energy are stored in the electric and gravitational fields 
and cancel each other out. The only energy is then that stored in the mass. 

We may think of this geometric energy as a "configuration energy" that arises from the spatial 
relationships of different parts of the spacetime to each other. By contrast in the next section 
where we include the gauge dependent terms we will see that our expression for the energy also 
includes "position" terms that arise due to the position of the different parts of the spacetime 
in the gauge potential. Of course, the form of the gauge potential is determined up to a gauge 
transformation by the matter so this view of the terms as being configurational versus positional 
is at best a rough way to think of them. 

4.1.2 Static total energy 

Next consider e + e m , the full energy density that was derived from the variational calculations. 
Then 

*- = / ^Me + n = =4±^ggg (82) 
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As we have emphasized before this expression is manifestly gauge dependent. Even worse however 
is the fact that this energy will in general diverge at the outer horizon of a black hole. Before we 
deal with that worry however, let us consider the usual r — * oo limit. 

If we demand that A a has the same spherical and time translation symmetries as the spacetime, 
then x = — &oot + f(r) where = lim r _ i . 00 $ is a constant and f(r) is an arbitrary function of 
r. Then, 

E to t-E = r- ^ ~l mr ~ E ° rd ^ ^(m + E^oo) + ^(m 2 + EO 2 + 2mE^ OCl ). (83) 
yV 2 - 2mr + El 2r 

Since the total energy is the sum of the geometric energy and the gauge dependent term, it isn't 
surprising that this Newtonian limit is the sum of the Newtonian limit of the geometric energy and 
the "positional" potential term. We can think of as the zero level of the potential throughout 
space (it remains even if Eq — > 0) and so by the thin shell analogy can think of -Eo^oo as the 
energy cost for creating charge in this spacetime before bringing it in from infinity. For extreme 
black holes recall that Eq £0 — E = m so the E tot — E_ = m + J* n d 2 x^e m and the only energy is 
the mass m plus the energy of the charge with respect to the potential. 

In most situations the exact choice of gauge is just a matter of convenience. For black hole 
spacetimes however, we noted above that for most gauge choices Etot will diverge on the horizon. 
This divergence can be directly traced to the fact that the Coulomb potential = — u a A a = 
-jpiSjj*- — dtx) also diverges at the horizon. To remove both divergences we must choose x such 

that dtx — > r as r — > r + where r + is the outer black hole horizon. That is we set the Coulomb 
potential to zero on the black hole horizon. Assuming that A a has the symmetries that we 
discussed above we must then choose = — Making that choice, after a little algebra we 
obtain 

It — 7"_i_ 

Etot = -rJ ± (84) 

V r — T- 

where r± = m ± \J m 2 — E 2 are the radial positions of inner and outer horizons. This gauge 
will also be used for the naked black holes. For extreme black holes r + = r_ = \Eq\ = m and 
so Etot = E everywhere. Physically we have chosen the gauge so that the potential energy is a 
constant and everywhere equal —m. The (in this case negative) electric potential energy cancels 
the mass-energy while at the same time the positive energy of the electric field cancels the negative 
binding energy of the gravitational field. 

Reflecting on this section we see that the total energy may in a very real sense be split into 
two parts. In the last section we saw that the geometric part depends only on the configuration 
of the spacetime. Examining the Newtonian limit we saw that it appears to include not only the 
gravitational but also the electromagnetic "configurational" energies. By contrast in this section 
we saw that the gauge dependent part exclusively deals with the potential of the matter relative 
to the gauge field. As we have seen, for a given solution to the Einstein-Maxwell equations the 
total QLE for a given surface may take any value (including zero) depending on the exact gauge 
choice. As such it is clear that this gauge dependent part of the energy is not reflected in the 
stress-energy tensor. On the other hand we should not then conclude that this gauge dependent 
part is entirely meaningless. It certainly plays a role equal to the geometric energy in both in 
thermodynamics M, M and in black hole pair creation [12|. 



4.1.3 Energies seen by radially moving observers 

We next consider the energies measured by spherically symmetric sets of observers who are moving 
radially in the RN spacetime. As before we take £lt to be r = constant, t = constant surfaces 
but now take T* a = N*u* a where u* a = \u a + r]n a = ^{u a + v\-h a ). As in section 3.4, v\- is the 
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speed of the T* a observers in the h a - 
evolving by u a = ~7f^?- 

Then, a straightforward calculation shows el = so 



Fd" direction as measured by a second set of observers 



E, 



Geo 



jE Geo = -jrVF. 



(85) 



Unfortunately, from the point of view of simplicity, within our gauge freedom e m I is not necessarily 
zero. Even if we restrict ourselves to gauge choices that give A a the same symmetries as the 
spacetime, we have x = ~&oot + f( r ) where as noted earlier $oo is a constant and / is an 
arbitrary function. Then 



47T 



Fd r f. 



(86) 



In the interests of simplicity however, we make the standard gauge choice for electrostatics and 
let d r f = 0. Then the Lorentz-type transformation laws will apply and 



TP* 



tot 



-ryf 



F 



Eq 
ry/F 



^ Eo 

$oc + — 

r 



(87) 



if we don't want this quantity to diverge at the horizon. 



As before we must choose 

To include the reference terms, we must first find a time vector T* a for the reference spacetime 
such that T* a T* = T* a T* and £t*£m8 = £T*&a/3- This prescription gives 
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1 - (1 - F)vlu a + v^VWhA , 



(88) 



where u a = df , h a = and t and r the usual time and radial coordinates for Minkowski space. 
Then ~ _ 

T* a n a v h VF 
Vy= ~ ~ ~ 



T* a u„ 



1 - (1 - F)v? 



and 7 = 71/ 1 — (1 — F)v, 



(89) 



Thus 



and we have 



E* = j_E = -nJl - (1 - F)vl, 
E* G eo -E*=r 7 ^1 - (1 - F)vl - VF^J , 



and 



E, 



tot 



E* =r 7 [Jl-{l-F)vt 



(Vf- e ° 




V ry/F 


r 



(90) 
(91) 

(92) 



As they stand these expressions are quite complicated and their physical interpretation isn't 
at all obvious. Thus let us consider the large r /small v limit. To first order in = and first order in 



771* 771* 

F Geo ~ F 



m 



2r 



m 



El) - -mvl, 



and 



E* tot -E* « (m + E $ c 



-(m + £ $ooK + 2~(^o + m2 + 2mE <Z> c 



(93) 
(94) 



These results are quite interesting. We see that the motion of the observers actually serves to 
decrease the quasilocal energy that they measure. Specifically we see that both the total and 
geometric boosted quasilocal energies are equal to their unboosted counterparts minus a kinetic 
term equal to the \ (Total Energy of Fields)^. This effect, also noted in [24|, is in stark contrast 
to both of the no-reference-term quantities whose measure increases with motion. Some discussion 
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of why this happens may be found in section 4.2.2 where we consider the equivalent effect for naked 
black holes. 

Next, in preparation for studying those naked black holes, we consider the specific case in which 
a spherically symmetric group of observers are falling into a RN black hole. These observers are 
assumed to have started with velocity close zero "close to infinity" and then fallen along a radial 
geodesic inwards towards the black hole. Rigorously the geodesic is the one that with respect to 
the standard time foliation has radial velocity zero at infinity and 1 at the outer horizon. Now, a 
test particle starting with velocity zero at radial coordinate ro and then allowed to fall towards a 
black hole on a radial geodesic will have coordinate velocity 



dr 

~d7 



-y/F(r )-F(r) 



(95) 



as a function of r, where r is the proper time. Thus an observer infalling on a geodesic that was 
static at infinity will have coordinate velocity ^ = — yl — F(r). 

1 „.a 



Let these observers measure time in the natural way (that is N = 1), T* 



Then the instantaneous radial velocity of the T* Q observers as measured in the static u a frame is 



and so the Lorentz factor is 7 = 



U/3 



(96) 



Substituting this value for 7 into equations (85,87,90) and making the gauge choice $00 = — ^ 



so that E\ ot doesn't diverge at the horizon, 

7-1* 

^Geo ~ 
7-1* 

E* -- 



-r, 



and 



-ry/2 - F. 



(97) 
(98) 
(99) 



Note that as r — > r + all of these take non-zero values. By contrast Eceo and Etot both are zero 
at r + . Also, keep in mind that for a near extreme black hole, r + ~ r_. Therefore for a black hole 
that is very close to being extreme, the observers will measure E\ ot to have a very large negative 
value as they approach the horizon. 
Including the reference terms, 



7-1* 771* 



and 



E, 



tot 



E* 



■(v^Tf-i), 



V2^F 



(100) 



(101) 



So near the horizon the infalling gravitational energy (including the reference term) goes to (V2- 
l)r_)_ compared to r + for the static gravitational energy. By contrast, the infalling total energy 
(including reference term) attains arbitrarily large negative values as the observers approach the 
horizon for black holes that are arbitrarily close to being extreme. Static observers however, will 
measure Etot — E = r + as they hover around the horizon. The difference is essentially due to 
the hugely boosted matter terms. As we have seen the boosting of the geometric terms has a 
comparatively minor effect. 
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4.1.4 The Hamiltonian 



Finally we return to static sets of observers measuring time according to the time coordinate t 
(that is lapse ./V = y/~F) and calculate the Hamiltonian energy measured by these observers. As 
usual the surfaces of observers are spherical t = constant and r = constant surfaces. Then 

H Geo = NE Geo = -rF, (102) 

H Geo -H = N(E Geo -E) = V^F(l-VF), (103) 

H tot = NE to t = -r + 2m + E dtX: and (104) 

Htat-H = N(E tot -E) = 2m + E d tX + V^F-r (105) 

771*^ | E^ TfiP" E^ 

In the large r limit, H Geo — H_ m m ^F -2 " an d H tot — H_m (m + Eb^oo) 2i~®~- Thus as is usual 

for asymptotically flat spacetimes the Hamiltonian corresponds to the QLE in the r — > oo limit. 
Note too though that in the large r limit these Hamiltonian's don't agree with the Newtonian 
limits that we discussed earlier. 

Finally let us illustrate the earlier comments on gauge invariance. To avoid the complications 
of singularities we redefine M to be the region of M contained by the two timelike hypersurfaces 
r = r\ and r = where r + < r\ < r%. Again we foliate that region according to the standard 
time coordinate t. Since we are only interested in the gauge invariance of the Hamiltonian and 
the reference terms are manifestly gauge invariant we ignore them for this calculation. 

Then the total Hamiltonian for the action IemcLU f° r the- region M is 

H M = XH tot = (n - r 2 ) + [dtxYrl (106) 
where = dtx\r=r 2 ~ &tx\r=ri- The sum is over the two boundary components. Consider the 



gauge dependence of this Hamiltonian. In section |3,2| we saw that we could only expect it to be 
gauge independent if M is was a region containing no singularities and £t^-ol = 0. Well, there are 
no singularities in M and a quick calculation shows that £r^a = implies that dtx is constant 
over E 4 . If this is true then [<9tx]£? = and the Hamiltonian is gauge independent as we expect. 

4.2 Naked black holes 

Finally, we consider naked black holes. They are low energy limit solutions to string theory and 
are characterized by the fact that static observers hovering close to their horizons see only very 
small curvatures while infalling observers are crushed by very large tidal forces. They are naked 
in the sense that even though they are not Planck scale themselves, Planck scale curvatures may 
still be observed outside of their horizons. Several classes of these holes were studied in a couple 



of papers by Horowitz and Ross [25| but here we will consider only those satisfying the equations 
of motion (^-^) . The naked black holes are then a subset of the following class of Maxwell-Dilaton 
black hole solutions. The metric is given by 

ds 2 = -F{r)dt 2 + -^-y + R{rf{d6 2 + sin 2 8d(p 2 ), (107) 

where 

F(r) = (r - r+ ^- r - ) and R(r) =r(l- , ( 108 ) 

and r + and r_ are respectively the inner and outer horizons of the black hole. The accompanying 
dilaton and electromagnetic fields are defined by 

m / r \2a/(i+a 2 ) 
e ~24, = ^ _ j ( 10Q ) 
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and 



G 

* F = —dt A dr. 



(110) 



These solutions are all magnetic black holes so as we have discussed earlier we work with the dual 
electromagnetic field tensor. The ADM mass and magnetic charge are 



r + 1 — a 2 r_ 



and 



Go 



r+r _ 



1/2 



(111) 
(112) 



Solving this pair of equations in terms of r + and r_ we find r± = j^i(M ± ^/M 2 - (1 - a 2 )G 2 ) 
for a 7^ 1 or r + = 2M and r_ = Gq/M for a = 1. Note that when a = these expressions reduce 
to those for a magnetically charged RN black hole. 

Massive near-extreme members of this class of solutions are naked. To see this we note that in 
terms of the orthonormal tetrad {u a , h a , § a , cp a } where u a = 1/VFdf, n a = \/Fd?, 6 a = 1/R9$ 
and (p a = l/(Rsin6)d® the non-zero components of the Riemann tensor are 
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TZ - - 


2' 




1-FR 2 


R 2 




= TZuif>u(p = 


TZ - - 


— T? — 



FR 

— ^, and 
2R ' 

FR FR 

~~2~R ~ ~R' 

In the above and the following, overdots indicate partial derivatives with respect to r. 



In the alternate moving tetrad {u a , n a , 9 a , (p a } where we recall that u a 



n a = jn 



r\n 



(113) 
(114) 
(115) 
(116) 

and 



rju 

components) 



the non-zero components of the Riemann tensor are (in terms of the non-moving 



TZ- - 



TZutfuif> 



TZ - - = — 



TZ 



T^wfiuip ~i~ f] (TZ U (p U (p -\- TZfiiphif)) 

2 



nipnip 



FR 

~m ~ 

FR 
~~2~R 



.FR 



and 



V 



,FR 
~R' 



(117) 
(118) 
(119) 



R 
FR 
~R 

T^0§00 is unchanged and we still have TZ-g-g = TZ^u^ and TZ n0ne = TZ ni p n $. Clearly if a = then 
R(r) = r and all of the components are the same as for the unboosted frame. 

Now if a and we define <5 = (1 — r_/r+) 1 ^ 1+a \ then the naked black holes are the subset 
of the above set whose parameters satisfy the conditions |^ *C -jjr ^ 1) where R+ = i?(r+). That 



is f » R 



which in turn is much larger than the Planck length. We note that if S <C 1 then 
r_ ~ r + and if R + 3> 1 then M, Go S> 1. Thus naked holes are near-extreme as well as being very 
large (relative to the Planck length). 



In the static frame as r — > r+, 



I TZuhuh | 



1 - 



TZ 



<p8ip8 
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Rj 



2r_ 



(l + a 2 )r+ 



« 1, 



(120) 
(121) 
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K U0U0 ( 1 - ^ , „, N _ ) < 1, and (122) 

\T^a<m\ -> tt^t ( 1 " m ,'1^- ) « L ( 123 ) 



1 / r_ 
2^ V 1 " (1+P)r + 

1 / r_ 
2^ V 1 ~ (T + "a 2 )r. 

Thus, all of the curvature components (and consequently the curvature invariants calculated from 
them) are small compared to the Planck scale. 

By contrast if we fix the tetrad to be that carried by the infalling observers considered in the 
previous example, rj 2 = j 2 v 2 = and as r — > r+, 

\n*m\ - (i + a 2)2 ^r^p2 » land (!24) 

a 2 r 2 1 

l^wl ~> ( 1 + a 2)2^^p2 >>L ( 125 ) 

Thus these infalling observers see Planck scale curvatures. Interpreting these components in terms 
of the relative acceleration of neighbouring geodesies we see that the observers are laterally crushed 
by huge tidal forces. 

4.2.1 QLE of naked black holes 

If we consider a spherical shell of observers falling into a naked black hole, we would expect these 
tidal forces to cause the area of the shell to shrink at a very rapid rate. Now rates of change of 
area factor largely in the formalism that we have developed in this paper. In particular e$ is, up 
to a normalization factor, exactly the (local) rate of change of the area of an infalling surface of 
observers. As such it is of interest to calculate the quasilocal energies measured by static versus 
infalling observers and to see how they compare to the observed curvatures. After a considerable 
amount of algebra we obtain the pleasantly simple 

- -ss^- < 127 > 

£ = -7^- (129) 

The gauge choice for the matter term is the same one that we used earlier. That is, we have chosen 
the gauge so that *A a \\ u a , as well as being static, spherically symmetric, and non-diverging on 
the black hole horizon. Though this is a long list of requirements, as noted earlier it is really 
little more than asserting that we make the standard gauge choice of electrostatics (or in this case 
magnetostatics) . 



For the infalling observers we again have T* a = ~^u a — y which implies v\- = — y/1 — F 

and 7 = l/s/F. By contrast the joint requirements that T*°T* = T* a T* and £T*g_ a p = £r*o'ai3 
give us 

T* a = y/ 1 + R 2 {1 - F)u a - Ry/T^Fn a , (130) 

Vy = , and (131) 

y/l + R 2 (l-F) 



7 = Jl + R 2 (l-F). (132) 
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Figure 2: Asymptotic and near horizon behaviour of the quasilocal energies for near-extreme dilaton-Maxwell black 
holes. We have 5 = (1 - r_/r + ) 1/(1+a2) < 1, i? + = r + 5 a2 > 1 and R\5 2 < 1, where i?+ = R(r + ). C\ is a constant on 
the same order as 1. 

Then, we have 



Eqeo = 


— y(r — r + ){r — r)R, 


(133) 


&Geo - 


— RR, 


(134) 


Etot = 


V r — r_ 


(135) 


TP* 

&tot - 


R 2 
r — r_ 


(136) 


E = 


—R, and 


(137) 


E* = 


-rJ\ + R 2 {1- F). 


(138) 



Evaluating these expressions at r = r + is straightforward with the only complication being R + = 
R(r + ). We have, 

*+ = tt^( 4 ° i + t)- < 139 > 

If a 2 <C J then the square of the coupling constant is extremely small even relative to 5, and 
R + ~ 1. In fact even if a 2 ~ 5 then i? + is of the same order as 1. By contrast for a 2 S> 5, 
+ ~ i+a 1 T ^ - Thus, we calculate the quasilocal energies for the cases a 2 <~ 5 (which 
includes the magnetic Reissner Nordstrom case for a = 0) and a 2 3> 5 separately. The results 
along with those for r — > oo are displayed in figure Note that for a 2 <~ 5 < 1, R + ~ r + . 

From figure § we see that static observers outside a naked black hole measure Ec e o, Etot —> 
near to the horizon while the infalling observers measure those same quantities to be very large. 
This effect occurs for both 5 <C a 2 and the a 2 <~ 5 (which include the RN holes) and so cannot be 
attributed to the "nakedness" of the holes. Of course since we have omitted the reference terms, 
both of these expressions blow up if we take the quasilocal surface out to infinity. 

If we include the reference terms, then near to the horizon Etot — E is very large for static 
observers where it is It is even larger in the absolute sense for infalling observers who 

measure it as —R+/5. Again however, those effects are seen by observers surrounding both naked 
and near-extreme RN holes and so cannot really be attributed to the curvature. As r — > oo the 
two expressions agree which is not surprising since as r — > oo the velocity of the infalling observers 
goes to zero. Note however that we do not obtain the ADM mass. 

More interesting are the measurements of Ea eo — E. If a ~ 1 and the holes are large {R\ S> 1) 
then we see that while static observers near to the horizon measure large values, sets of observers 
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falling into naked black holes actually measure very small values for this quasilocal energy. By 
contrast observers falling into an RN hole will measure large values. In fact we can see that for 
a ~ 1 and 3> 1 then these infalling observers will measure Eq — E* <C 1, if and only if the 
black hole is naked. Thus this is an alternate characterizing feature of naked black holes when 
the coupling constant is of a reasonable size. The equivalence is broken if a 2 <~ 5 in which 
case the static and infalling observers both measure large energies. Consider for example the case 
where a 2 = 5. Then the black hole can still be naked if S (and therefore a 2 ) is small enough that 
5 2 R 2 « 1. 



4.2.2 Why do naked holes behave this way? 

At the beginning of the previous section we suggested that the curvature results could be under- 
stood in terms of the rates of change of the surface area of shells of infalling observers. In this 
section we explore this idea if more detail and also use it to provide an explanation of the Ec eo — E 
result. 

First we quantify the expectation that the surface area of a shell of infalling observers will be 
changing extremely quickly as they cross the horizon of a naked black hole. Recall that naked 
black holes are near extreme and so the singularity sits "just behind" the horizon (r_ ~ r+). More 
rigorously, Horowitz and Ross p^j noted that an observer passing through the horizon after falling 
from ro (the situation described by equation (]95[)), will hit the singularity at r_ after a proper 
time of As <~ r +^^~ = ■ Thus a set of observers infalling on geodesies that were stationary 
at infinity (-F(ro) = 1) will only have a very short time before they reach r_. At r_, R(r) — > and 
so the area of the shell goes to zero. However, by assumption R + ^ 1 so at the horizon itself, that 
same area is very large. For the area to go from very large to zero in such a small time, we would 
naively expect it to be decreasing very quickly as the observers pass the horizon. To quantify this 



we use (56) to show that the fractional rate of change of the area of the surface Qt as measured 



by the observers who inhabit that surface is 



A 8vr L d 2 x^fde^ 2R+ 2 ( „„ 2 a 2 



A j Qt d 2 x^ R + (1 + a 2 )R + ^ + 5 2 ) ' (M0) 

If a 2 3> 5 (that is it isn't pathologically small), ^ ~ > 1 as we would expect. By contrast 

for the RN case (a = 0), ~ -C 1 However, our expectations are confounded if 5 >~ a 2 ^ 
in which case the hole remains naked even while the rate of change is more along the lines of the 
RN values. In that case the extremely small value of a suppresses the rapid decrease in area until 
the observers get even closer to the singularity (basically r — r_ <C a 2 ). 

These rates of change of the area also nicely explain why Eo eo — E is small while the observed 
curvature components are large. Recall that to define the reference term E we had to embed 
into flat space along with a vector field T_ a defined so that if we evolved Of with that vector field 
and made only intrinsic observations on the resulting timelike three surface then observers could 
not tell whether they were in the original or reference spacetimes. In particular, the area of 
should change at the same rate. Thus, if the area decreases extremely rapidly, the embedded shell 
of observers in the reference spacetime would have to be moving at a correspondingly fast speed. 
Equation ( |131| ) quantifies this telling us that Vy = R/y 1 + R 2 at the horizon. Then for a 2 ~ 1, 
R 3> 1 =>- Uj_ ~ 1 and the observers have to move at close to the speed of light in the reference 
time to match the rate of change of the area. By contrast, for a 2 ~ 0, R ~ 1 =>• Vy ~ |. The 
area is changing at a relatively leisurely rate so the observers do not need to move so fast in the 
reference time. 

For observers moving at extremely rapid velocities there is a sense in which the relativistic 
effects of their speed become more important than those due to gravity. To see this recall equation 
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(3.4) in which we saw that e 2 — el 2 is a constant independent of the speed of the observers. Now, 
by construction e"l is the same in both the reference and original spacetime and so we can rewrite, 

Eheo ~ K* = d 2 x^ (yJg+iV 2 - Ve 2 + el 2 ) • (141) 
If el is much larger than e and e, then at the horizon 

Eheo ~ M* » ~ [ d 2 xV^ (^A) , (142) 



2 Jn t V el , 

and so we see that as el becomes larger and large the observers quasilocal energy becomes smaller 
and smaller. Physically, though e* and e* are boosted to be very large, the difference between 
them simultaneously becomes smaller and smaller. In particular for naked black holes we have 

Eheo ~E*~ 2nR 2 + ( pi = = — (143) 



V 2R+ " A' 

and we see that in this case the geometric quasilocal energy is actually the inverse of the rate of 
change of the area. 

By contrast E^ ot — E_* includes matter terms which are also boosted to be very large. There 
is no corresponding term in the reference spacetime to cancel these large terms out. The result is 
the matter terms dominate over the geometrical terms in E^ ot — E_* and so this total quasilocal 
energy is very large. 



5 Discussion 

In order to investigate how different sets of observers would define energy, angular momentum 
and charge in systems of finite size, it is necessary to make use of a formalism which explicates 
the relationships between those observers. The non-orthogonal formalism developed in previous 

is designed to address this issue. We have in this paper extended it to 
include electromagnetic and dilatonic matter fields. 

Our applications of this extended formalism have yielded a number of results. The QLE 
naturally breaks up into a gauge independent geometric part and a gauge dependent part. Likewise, 
the Hamiltonian only inherits a subset of the gauge invariance of the full action. These situations 
arise because of the manner in which the quasilocal formalism separates out the bulk and boundary 
terms of the action. For a particular gauge choice which is natural to make in stationary spacetimes, 
the gauge dependence of the Hamiltonian is eliminated and that of the QLE is reduced to that 
resulting from where we choose to set the zero of a the Coulomb potential. 

We have also shown that the quasilocal energy, angular momentum, and spatial stress have the 
expected covariance properties under Lorentz-type transformations, allowing us to meaningfully 
relate observations made by static observers to those made by non-static ones. However this 
covariance is not respected by the reference background spacetimes that are used to render the 
various quasilocal quantities finite. The reason that the covariance is not respected is that the 
quasilocal quantities in the reference spacetime transform with respect to a different velocity than 
the velocity of the quasilocal surface in the original spacetime. The requirement that the two 
surface evolve in the same way in both spacetimes (which after all have different geometries) 
forces these two velocities to be different. 

When we apply this formalism to various physical situations, we find a variety of results which 
both confirm and confound our intuition. The thin-shell formalism provides a natural operational 
definition of the QLE: the QLE contained within a 2-surface can be defined as the energy of a 
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shell of matter that has the same intrinsic geometry as this 2-surface and a surface stress energy 
tensor defined so that the spacetime outside the 2-surface is identical to the original spacetime, 
while inside it is identical to the reference spacetime. A consideration of RN and naked black holes 
yields an interesting set of similarities and differences between the QLEs in each case as measured 
by differing sets of observers. For both types of black holes, the geometric energy EQ eo and 
total energy Etot vanish for static observers near the event horizon, whereas infalling observers 
measure these same quantities to be very large. When the reference terms are included, both 
static and infalling observers measure Etot — E to be very large for both kinds of holes. However 
measurements of EQ eo — E_ will differentiate between the two kinds of holes. Infalling observers will 
find this quantity to be large for RN holes but small for naked holes, provided that the coupling 
constant a is not too small. 

Some outstanding issues remain. The inclusion of dyons remains an unsolved problem within 
our formalism. A natural first approach for dealing with such objects would be to consider different 
patches of the quasilocal surface, with different (but gauge-related) nonsingular gauge fields defined 
on each. Unfortunately this is not sufficient to solve the problem. A simple way of seeing this is 
to note that for a pure Maxwell spacetime F 2 = 2(B 2 — E 2 ). With E a and B a not being treated 
on an equal footing in the action it is not so surprising that it they are not treated on an equal 
footing in the QLE defined from that action. Thus it appears that the action itself will have to 
be modified in some way to resolve this problem. 

The relationship between the thin-shell formalism and QLE defined with AdS/CFT-inspired 
intrinsic reference terms also deserves further study. Although the correspondence between the 
two formalisms does not carry through, the intrinsic reference terms used so far have only been 
those which are the least divergent for quasilocal surfaces of large mean radius. Other intrinsic 
reference terms could be included for finite-sized quasilocal surfaces - their role in the thin-shell 
QLE correspondence remains to be explored. 
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